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ABSTRACT 

The  equations  of  motion  of  a  puck  sliding  frictionlessly  inside  a 
parabolic  bowl  can  be  straightforwardly  deduced  using  the 
conservation  laws  of  mechanical  energy  and  angular  momentum.  But 
the  solution  of  these  equations  requires  that  they  be  recast  into  the  form 
of  Newton’s  second  law.  The  simple  example  of  a  ball  in  vertical 
freefall  illustrates  why  this  is  necessary  and  how  to  perform  the 
conversion.  The  method  is  then  applied  to  the  richer  problem  of  a  puck 
gliding  on  a  paraboloidal  surface  for  which  the  nonlinear  equations 
require  numerical  solution.  A  rich  variety  of  orbital  patterns  of  the  puck 
is  found. 


Introductory  Example  of  One-Dimensional  Freefall 

Consider  a  ball  thrown  straight  upward  (which  will  be  designated  as 
the  +z  direction)  from  the  origin  with  an  initial  velocity  of  v0z.  Let’s  find 
its  resulting  path  of  motion  z(t)  in  the  absence  of  air  resistance.  Because 
mechanical  energy  is  conserved  (for  the  system  of  ball  and  earth),  the  sum 
of  the  kinetic  ( K)  and  gravitational  potential  (U)  energies  at  any  point  on 
the  ball’s  path  can  be  written  as 

k  +  u  =  kq  +  u0,  (1) 

where  the  subscript  “0”  throughout  this  article  denotes  the  initial  instant 
t  =  0  .  Choosing  the  gravitational  reference  position  to  be  at  the  origin  and 
assuming  the  ball’s  altitude  never  gets  large  compared  to  Earth’s  radius, 
Eq.  (1)  becomes 
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12  12 
-  mv  +  mgz  =  -  mu, 

2  z  6  2  Os 


where  m  is  the  mass  of  the  ball,  g  =  9.80  N/kg  is  Earth’s  surface 
gravitational  field,  and  v_  =dz/  dt  is  the  velocity  of  the  ball.  Equation  (2) 
can  be  rearranged  as 


2  gz. 


(3) 


Unfortunately  this  equation  is  double-valued  and  cannot  be  uniquely 
solved  as  written.  At  any  given  height  z,  there  are  two  solutions,  one 
corresponding  to  the  ball  traveling  upward  with  a  positive  velocity  and  the 
other  to  the  ball  descending  with  an  equal-magnitude  negative  velocity.  In 
order  to  circumvent  this  ambiguity,  the  time  derivative  of  Eq.  (3)  can  be 
taken  to  produce  the  readily  solvable  form 


(  dz)  d  z  \  _  f  dz 

2  -r  -^r  =-2  s  — 


a  ——2 

z  ^ 


2  2 

where  a_  =  d  z  /  dt  is  the  acceleration  of  the  ball.  The  final  equation  is 
simply  Newton’s  second  law  with  the  ball’s  mass  divided  out  of  both 
sides.  Integrating  it  twice  with  respect  to  time  gives  the  expected  solution 
z  =  vQzt-±gt2. 

In  this  easy  example,  one  could  alternatively  solve  Eq.  (3)  by  manually 
changing  the  sign  of  the  square  root  of  the  right-hand  side  of  the  equation 
after  the  topmost  point  of  the  trajectory  is  reached  by  the  ball.  But  this 
procedure  becomes  cumbersome  if  the  orbit  has  a  large  number  of  turning 
points.  In  such  a  case,  it  is  easier  to  differentiate  the  energy  equation  with 
respect  to  time  and  then  solve  the  resulting  second-order  equation,  as  was 
done  above.1  Let’s  now  apply  this  method  to  the  richer  problem  of  interest 
in  this  paper. 

Orbiting  On  a  Frictionless  Parabolic  Surface 

Suppose  that  a  puck  is  sliding  frictionlessly  about  the  bottom  of  a 
concave  bowl  which  has  cylindrical  symmetry  around  the  vertical  axis  z, 
described  by  the  parabolic  cross-sectional  profile 
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z  =  l2kp2  (5) 

using  cylindrical  coordinates,  p,  (j),  z,  as  illustrated  in  Fig.  1.  The  origin  of 
the  coordinate  system  is  at  the  vertex  of  the  bowl,  and  a  factor  of  lA  has 
been  included  in  Eq.  (5)  to  avoid  factors  of  2  that  otherwise  arise. 


Fig.  1.  Free-body  diagram  indicating  the  normal  (N)  and  gravitational 
forces  (mg)  acting  on  the  puck  (indicated  by  the  dot)  when  it  is  located 
at  arbitrary  position  ( p,(j),z )  •  The  paraboloidal  surface  has  slope  tan# 
in  the  radial  direction. 


Energy  conservation  implies  that 

\mv2  +  mgz  =  constant  =>  v2  +  gkp2  =  constant  (6) 

where  v  =  v~  +  V2  +  v  and  the  first  constant  has  been  divided  by  a 
factor  of  Vi  m  to  get  the  second  constant.  Here  vp  =  clp  /  clt ,  =  pco 

(where  co  =  d(f)/dt  is  the  puck’s  angular  velocity  about  the  axis  of 
symmetry),  and  v_  =  dz  /  dt  =  kp  dp  /  dt .  Since  neither  gravity  nor  the 
normal  force  exerts  a  vertical  torque  on  the  puck  about  the  origin,  the  z- 
component  of  the  angular  momentum  is  constant  and  therefore  equals  its 
initial  value, 


L 

z 


2  2 

=  mp  (O  =  mp0a>0 


CO  = 


(  \2 

Po 


V 


CO, 


0- 


J 


(7) 


Inserting  this  expression  into  the  speed  squared  in  Eq.  (6)  and  taking  the 
time  derivative  to  eliminate  the  constant  yields 
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d_ 

dt 


(l  +  k2p2^v2  +  PqCOqP  2  +  gkp2 


0. 


(8) 


The  derivative  can  be  performed  and  a  factor  of  2v  divided  out  of  every 
term,  in  analogy  to  how  Eq.  (4)  was  obtained  from  Eq.  (3),  to  get 


4  2-3 

p0co-p 


kp^g + kv 


1  +  k2p2 


(9) 


where  a  =  d2p  /  dt 2 .  This  equation  can  also  be  obtained  (but  with 
considerably  more  effort)  by  finding  the  two  orthogonal  surface  tangential 
components  (to  avoid  the  unknown  normal  force)  of  Newton’s  second  law 
in  cylindrical  coordinates. 

One  final  step  is  helpful  before  proceeding  to  a  computer  solution. 
Equation  (9)  can  be  rewritten  in  terms  of  the  dimensionless  variables 
R  =  kp  and  T  =  coQt  as 


d2R 

dT2 


C  +  (dR  /  dT f 


(10) 


2 

where  C  =  gk  /  C0q  is  a  dimensionless  constant.  This  is  a  second-order 
differential  equation  to  be  solved  with  the  initial  conditions 
R( 0)  =  Rq  =  kpQ  and  V(0)  =  V0  =  kvpQ  /  coQ  where  V  =  dR  /  dT .  Suppose 
the  initial  angular  velocity  is  chosen  so  that  the  puck  travels  in  a  stable 
counter-clockwise  circular  orbit  around  the  vertex  of  the  bowl.  The  puck  is 
then  given  a  quick  push  toward  the  rim  of  the  dish.  The  push  provides  a 
radial  impulse  to  the  puck.  (Note  that  a  radial  impulse  does  not  change  the 
value  of  Lz.)  Prior  to  the  push,  R  must  have  the  constant  value  R0  so  that 
dR  /  dT  and  d  R/  dT  are  both  zero,  and  Eq.  (10)  therefore  implies  that 
C  =  1 .  In  turn  this  result  requires  that  coQ  =  ( gk )  regardless  of  the 
puck’s  position  on  the  surface.  This  is  a  special  property  of  a  parabolic 
dish  and  is  the  reason  that  the  surface  of  a  rotating  liquid  settles  into  a 
paraboloidal  shape,  a  property  that  can  be  exploited  to  make  the  primary 
collecting  mirror  of  a  reflecting  telescope.2 

Once  Eq.  (10)  is  solved  for  R(T ),  it  can  be  substituted  into  Eq.  (7) 
written  in  the  dimensionless  form  dp  /  dT  =  ( R{)  /  R)  .  That  result  can 
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then  be  integrated  to  obtain  <f)(T )  with  the  initial  condition  0(0)  =  0  (by 
choosing  the  x-axis  to  point  to  the  puck’s  position  at  the  instant  of 
application  of  the  radial  impulse).  The  results  can  then  be  plotted 
parametrically  to  give  an  overhead  view  of  the  xy-coordinates  of  the  puck 
in  the  dimensionless  form 

X  =  Rc  os0  and  F  =  /?sin0.  (11) 

Here  is  the  complete  code  we  wrote  to  solve  and  plot  the  motion  of  the 
puck  using  the  commercial  software  program  Maple™  for  the  case 
Rq  =  1  =  VQ ,  as  graphed  in  Fig.  2(a): 

R0:=1;  V0:=1; 

eqR:=diff(R(T),T,T)=(R0A4-R(T)A4*(l+diff(R(T),T)A2))/R(T)A3/(l+R(T)A2); 

eqphi:=diff(phi(T),T)=(R0/R(T))A2; 

sol:=dsolve(  { eqR,eqphi,R(0)=R0,phi(0)=0,D(R)(0)=V0 } ,  { R(T),phi(T) }  .numeric); 

r:=T->rhs(sol(T)[2]);p:=T->rhs(sol(T)[4]); 

X:=T->r(T)*cos(p(T));  Y :=T->r(T)*sin(p(T)); 

plot(['X(T),,,Y(T),,T=0..50*Pi],scaling=constrained); 

By  varying  the  initial  values  R0  and  V0  in  the  first  line,  a  rich  variety  of 
orbital  patterns  result;  two  further  examples  are  plotted  in  panels  (b)  and 
(c)  of  Fig.  2,  chosen  to  illustrate  some  common  patterns.  Our  school  has  a 
site  license  for  Maple™  and  students  are  introduced  to  its  use  in  their 
introductory  calculus  sequence  and  could  be  given  the  above  code  with 
which  to  experiment.  At  other  schools,  Mathematica™  or  implementation 
of  Euler’s  method  in  a  spreadsheet  such  as  Excel™  might  be  a  better 
choice.3  However  the  comparative  simplicity  of  the  code  above  makes  this 
a  good  example  with  which  to  introduce  students  to  algorithmic  software 
packages. 

Further  insight  into  the  puck’ s  motion  is  obtained  by  making  the  radial 
impulse  very  weak,  so  that  the  circular  orbit  is  only  slightly  perturbed.4  In 
that  case,  it  is  easier  to  see  the  resulting  small  effect  by  jumping  into  a 
frame  of  reference  that  rotates  with  the  puck’s  initial  angular  speed  of  co0. 
The  xy-coordinates  of  the  puck  in  this  rotating  frame  can  be  computed 
using  Eq.  (11)  provided  we  replace  0  by  0  -  (0Qt  =  0  -  T  .  An  example  is 
plotted  in  Fig.  2(d).  The  puck  starts  on  the  x-axis  at  (i?0,0)  and  travels5 
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clockwise  with  very  nearly  uniform  circular  motion  of  dimensionless 
diameter  V0  at  an  angular  frequency  of  2 co0.  That  is,  the  puck  performs 
one  clockwise  orbit  in  the  rotating  frame  during  the  time  that  the  puck 
rotates  counter-clockwise  halfway  around  the  bowl  in  the  lab  frame.  This 
trajectory  is  a  result  of  the  Coriolis  force  which  produces  a  rightward 
deflection  of  the  puck  in  the  rotating  frame,6  analogous  to  the  rotation  of 
hurricanes  in  the  northern  hemisphere  of  the  earth.  The  radially  outward 
centrifugal  force  is  almost  perfectly  canceled  by  the  inward  component  of 
the  normal  force. 


Fig.  2.  Overhead  views  of  the  trajectory  of  the  puck  (a)  in  the  lab  frame 
for  Rq  =  1  and  Vq  =  1  over  the  interval  0  <  T  <  50 n  ;  (b)  in  the  lab 
frame  for  Rf)  =  1  and  Vq  =  8  over  the  interval  0  <  T  <  1 5()7T  ;  (c)  in  the 
lab  frame  for  R0  =  0.05  and  VQ  =  0.5  over  the  interval  0  <  T  <  25k  ; 
(d)  in  the  rotating  frame  for  Rq  =  0.01  and  V()  =  0.0001  over  the 
interval  0  <  T  <n  (in  a  highly  magnified  view). 
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For  Further  Investigation 

At  least  two  interesting  lines  of  inquiry  are  left  for  future  work,  the 
first  theoretical  and  the  second  experimental: 

(1)  Under  what  circumstances  is  the  orbital  motion  closed?  Close 
inspection  of  Fig.  2(b)  indicates  that  the  orbit  appears  to  repeat  after 
tracing  out  19  lobes.  In  contrast,  the  pattern  in  Fig.  2(a)  is  starting  over 
(after  25  time  periods  of  2n  /  coQ )  at  a  slightly  shifted  angular  position.  By 
writing  V  =  clR  /  dT  =  ( R()  /  R)  ( clR  /  d(p)  and  equating  it  to  the  positive 
square  root  of  V  from  Eq.  (6)  as  the  puck  travels  from  closest  to  farthest 
approach  from  the  bowl’s  vertex,  one  can  integrate  to  find  an  expression 
for  A0  along  this  path.  The  orbit  is  closed  if  n  /  A(f)  is  a  rational  number. 
(In  particular  if  that  number  is  an  integer,  then  the  orbit  never  crosses 
itself.)  Similarly,  Eq.  (10)  can  be  recast  into  an  orbital  differential  equation 
for  R((f>)  rather  than  R(T ) . 

(2)  To  investigate  experimentally  the  trajectories  described  in  this  paper, 
one  could  construct  a  parabolic  “air  hockey”  table  by  drilling  holes  in  a 
suitable  dish  and  blowing  air  through  them.  Alternatively  one  could  roll  a 
marble  on  an  old  parabolic  mirror  or  satellite  television  dish  and  modify 
the  present  theory  to  include  frictional  forces.  (One  could  even  spin  the 
dish  to  keep  the  marble  from  slowing  down.)  For  comparison,  interesting 
effects  occur  when  a  ball  rolls  without  slipping  on  the  surface  of  a  rotating 
flat  plate,7  on  the  inner  surface  of  a  vertical  cylinder  such  as  a  golf  cup,8 
on  the  surface  of  an  elastic  membrane,9  or  on  the  inner  surface  of  a 
sphere.10 
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